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Abstract 



cn A methodology to calculate the friction coefficient of an aggregate in the 

,— I continuum regime is proposed. The friction coefficient and the monomer 

^ shielding factors, aggregate- average or individual, are related to the molecule- 

aggregate collision rate that is obtained from the molecular diffusion equa- 
tion with an absorbing boundary condition on the aggregate surface. Cal- 
culated friction coefficients of straight chains are in very good agreement 
with previous results, suggesting that the friction coefficients may be ac- 
curately calculated from the product of the collision rate and an average 



O momentum transfer, the latter being independent of aggregate morphology. 

I ^1 Langevin-dynamics simulations show that the diffusive motion of straight- 

chain aggregates may be described either by a monomer-dependent or an 
^ aggregate- average random force, if the shielding factors are appropriately 

^—1 chosen. 

O 

00 Keywords: Straight-chain aggregates, friction coefficient, diffusion 

^ coefficient, collision rate, Langevin simulations 

1^ 

O 

, I 1. Introduction 

.1^ The hydrodynamic drag on a fractal aggregate suspended in a viscous 

rS fluid determines many of its dynamical properties, such as sedimentation, 

d agglomeration, and mobility (diffusive, mechanical, electrical) [1]. Fractal 
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aggregates arise from the agglomeration of smaller, primary spherules, here- 
after called monomers, that do not coalesce, but rather retain their identity in 
the resulting aggregate. The calculation of the drag force can be a formidable 
task, often requiring simplifying assumptions like spherical [21 [3] or ellipsoidal 
symmetry |3]. 

The friction coefficient of non-spherical aggregates in the continuum regime 
(fluid mean free path much smaller than the monomer radius Ri) has been 
calculated by treating the fractal aggregate as a porous particle in a viscous 
creeping flow (Stokes flow) of constant [1] or variable [21 |3l |5] permeability. 
Alternatively, Filippov [B] performed a multipole expansion of the Stokes flow 
velocity with slip and non-slip boundary conditions on the aggregate surface. 
Happel and Brenner [7] calculated the hydrodynamic forces on a collection 
of spheres by the method of reflections obtaining a slowly converging series. 
In these calculations the hydrodynamic force was obtained by integrating 
the fluid stress tensor over the aggregate surface. For straight chains, Dah- 
neke [S], based on an extrapolation of experimental data, proposed an empir- 
ical correlation for the friction forces parallel and perpendicular to the chain 
symmetry axis. 

Herein, we propose a versatile, albeit approximate, procedure to deter- 
mine the friction coefficient of an aggregate in creeping flows in the continuum 
regime. We argue that the friction coefficient of a fractal-like aggregate may 
be related to the ratio of two molecular collision rates: the molecule-aggregate 
collision rate and the molecule-monomer collision rate. These molecular col- 
lision rates are calculated from the diffusion equation for the fluid density 
with appropriate boundary conditions. Hence, we suggest that the friction 
coefficient of fractal-like objects may be accurately (but approximately) cal- 
culated from the solution of the Laplace equation (with appropriate boundary 
conditions) for the fluid density, without requiring the solution of the Stokes 
equations for the fluid velocity. 

The relationship between collision rates and friction coefficient allows us 
to introduce the monomer shielding factor (either individual or aggregate- 
average) that provides a measure of monomer shielding within an aggregate. 
We apply the proposed methodology to calculate the friction coefficients of 
straight-chain aggregates. The calculated friction coefficients are compared 
to previous analytical results and numerical results to justify the approxi- 
mations made a posteriori, and to validate the proposed methodology. We 
remark, however, that the methodology is general enough to be applicable 
to the calculation of friction coefficient of a general fractal-like aggregate. 



Calculated shielding factors of individual monomers in straight chains 
composed of fc = 5, 8 monomers are used to obtain the chain diffusion co- 
efficient by Langevin-dynamics simulations. These coefficients are found to 
be in very good agreement with the diffusion coefficients obtained from the 
ratio of molecular collision rates. 

Such Langevin simulations of aggregate formation and motion, which de- 
pend explicitly on the inter-monomer interaction potential, offer a potentially 
useful tool to investigate the effect of inter-monomer forces on agglomeration 
dynamics and the shape of the resulting aggregate. Previous works have sug- 
gested, either experimentally or numerically, the importance of inter-particle 
forces. For example, Chakrabarty el al. [H] hypothesized that electic forces 
were responsible for the unexpectedly low fractal dimension (between 1.2 and 
1.5) of a minority of soot aggregates formed in a premixed ethene-oxygen 
flame, di Stasio et al. [1^, in their experimental and numerical studies of the 
kinetics of agglomeration and growth of soot nanoparticles in an ethylene-air 
diffusion flame, suggested that the number of collisions is enhanced due to 
the attractive van der Waals forces. Kostoglou and co-workers [HI [12] stud- 
ied theoretically the time-dependent evolution of aggregate morphology, as 
described by the aggregate fractal dimension. Their analyses was based on 
a proposed constitutive law for the agglomeration of two aggregates, a law 
that depends on the inter-monomer (and hence inter-aggregate) interaction 
potential. The Langevin simulations described herein and in, e.g., Ref. [13], 
could address these theoretical and experimental observations by providing 
a technique to validate them. 

2. Molecule-aggregate collision rate and friction coefficient 

2.1. Methodology 

The motivation of our work stems from experimental measurements of 
the mass transfer (attachment) coefficient of gas molecules to nanoparticles 
of a variety of species, shapes, and sizes [lU [15]. The rate of transfer of 
gas molecules of density p to identical aggregates composed of k monomers 
(fc-aggregates) of concentration A^^ is 

f^=-K,N,p , (1) 

where K^. is the total mass transfer coefficient of fluid molecules per aggre- 
gate, namely Kk = hkSk where hk is the average heat transfer coefficient 



as usually defined in the chemical engineering literature (per unit surface 
area) and Sk the geometric surface area of the aggregate. Equation ([I]) is 
valid for a molecule-aggregate sticking coefficient of unity. Such mass trans- 
fer measurements may be performed by labelling atoms or molecules, via, for 
example, radioactive labelling (epiphaniometer) or electrical charging (diffu- 
sion charger), and subsequently detecting the atoms or molecules attached 
to the aggregate. Siegmann and Siegmann |Tlj, who measured total mass 
transfer coefficients over a wide range of aggregate mobilities via diffusive 
electrical charging, argued that the product of the mass transfer coefficient 
times the aggregate (electrical) mobihty B^ is approximately constant 

Kk X Bic = constant , (2) 

where Bj., the aggregate mobility, is the ratio of aggregate terminal velocity 
to the steady-state force causing it [TB]. 

Keller et al. [15] presented further experimental measurements support- 
ing the claim that Eq. ^ hold (within experimental error) for a variety of 
experimental conditions and aggregate shapes. They suggested that Eq. ^ 
constitutes an "empirical scaling law". They justified it by arguing that 
both the mass transfer coefficient and the aggregate mobility depend on the 
aggregate surface area exposed to the fluid. They referred to the particle sur- 
face area upon which mass, momentum, and energy transfer are dominant 
as their active surface area. Accordingly, the mass transfer coefficient is pro- 
portional to the particle active surface area, whereas the mobility is inversely 
proportional rendering their product a constant [HI [T3] . In particular, the 
total mass transfer coefficient K^. to a sphere in the continuum regime is 
proportional to its diameter dp, whereas in the free molecular regime it is 
proportional to its geometric surface area (~ dp). 

If multiple scattering events are neglected (sticking coefficient of unity), 
the attachment rate K^Nk is proportional to the molecule-aggregate colli- 
sion rate. For Stokes drag the mobility of a fc-aggregate is inversely pro- 
portional [IE] to the aggregate friction coefficient fk- Equation (|2]), then, 
implies that the collision rate Kk between a molecule and a fc-aggregate 
is proportional to the friction coefficient, Kk = cfk- The proportionality 
constant may be evaluated by considering, as experimental measurements 
suggest, that Eq. ^ is approximately valid for fc = 1. At steady state and 
for a spherical monomer, the integral molecule-monomer collision rate Ki is 
calculated from the the gas diffusion equation with an absorbing boundary 



condition on the monomer surface. Specifically, 



Ki= 3i-dS= Ji-sdS = ATTDgRip^ , (3) 

Jsi Jsi 

where the diffusive flux is Ji = —Dg Vp with p{r) = poo(l — Ri/r), poo is the 
gas-phase (fluid) density far away from the aggregate, and Dg the gas self- 
diffusion coefficient (see, for example, Ref. [1]). The monomer surface area is 
denoted by S*!, s is the unit vector orthogonal to the monomer surface, and 
dS the surface element [dS = R\ sin 9d9 d(f)). 

The Stokes friction coefficient of an isolated monomer is /i = rriiPi = 
GirpgRi with pg the fluid viscosity, leading to the proportionality constant 
c = 2DgPoo/{3pg) (for a similar calculation of the particle active surface area 
in the continuum regime cf. Ref. [E]). It is more convenient to eliminate 
the proportionality constant by considering the drag force fk relative to the 
total drag on k isolated monomers /i, 

^ = A f4) 

kK, kf\ - ^ ^ 

Theoretical arguments provide partial support of Eqs. (|2]), Q. In kinetic 
theory the mean force exerted on a particle by fluid molecules is the product 
of the molecule-particle collision rate times the mean momentum transfer per 
collision [TH]. The main approximation of our work, and the approximation 
implicit in Eqs. (pi), (|4]), is that the dominant contribution to the friction force 
arises from the collision rate. If the mean momentum transfer is taken to 
be independent of aggregate morphology the ratio of the friction coefficients 
fk/{kfi) becomes the appropriate ratio of collision rates shown in Eq. dll). 
It is, thus, apparent that Eq. dil) is approximate. The friction force arises 
from the momentum transfer during molecule-particle collisions; even if the 
sticking probability is unity not all collisions transfer the same momentum 
since gas velocities are distributed according to a probability density function 
that may, in principle, be calculated from the Boltzmann equation. 

The Stokes friction coefficient of a /c-aggregate is expressed as fk = krriiPk 
where Pk is the average friction coefficient per unit monomer mass and kirti 
the aggregate mass [13]. Equation ^ then leads to our main result that 
relates the friction coefficient of a general-shaped aggregate to the ratio of 
two molecular collision rates 

kR-rj.^'"' ■ ^'^ 



The last equality defines the average monomer shielding factor 77^ as the 
ratio of the average friction coefficient of a monomer in an aggregate to the 
friction coefficient of an isolated monomer. The shielding factor provides a 
measure of the shielding of a monomer by other monomers in an aggregate; 
as such, it depends on aggregate morphology. It has been used to calculate 
corrections of the Stokes drag on a porous medium consisting of identical 
non-interacting spheres [I9] , in modifications of heat transfer to an aggregate 
due to monomer shielding [20], and it has been related to the aggregate 
diffusion coefficient [13]. The friction coefficient, calculated according to 
Eq. (|5|, naturally leads to the Stokes-Einstein diffusion coefficient D^ of a 
A;-aggregate via 

kmiPk ^ krik 

where the Stokes-Einstein monomer diffusion coefficient is -Di[= kBT / {mifii)]. 

Furthermore, the aggregate mobility radius, defined as D^ = ksT /{Qn^gRk), 

is 

R, 

k rjk . (7) 



Rk 



Ri 

For the ideal aggregates defined in Ref. [13], aggregates in the free-draining 
approximation whereby the hydrodynamic forces on a monomer are inde- 
pendent of its state of aggregation, 77^ = 1. The shielding factor also gives 
the ratio of a fc-aggregate active surface area to the active surface area of k 
isolated monomers. 

The dynamic shape factor of particles of arbitrary shape, a correction 
factor used to account for the effect of shape on particle motion, is also 
related to the average monomer shielding factor. The aggregate dynamic 
shape factor Xs is defined by 

^^ = a ^ ' (^) 

bTT^gReq 

where the Req is the radius of the equivalent volume sphere, Req = k^^^Ri. 
Hence, 

Xs = %fc2/3 . (9) 

As in the case of an isolated monomer, the integral collision rate may 
be calculated from the fluid diffusion equation in steady state [V^p(r) = 0], 
with an absorbing boundary condition on the aggregate surface [p(rsur) = 0, 
neglect of multiple scattering events] and constant density far away from the 
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aggregate, (p -^ poo for |r| — )■ oo). The diffusive ffux J^ = —DgVp yields 
the coUision rate for a generic /c-aggregate 



Kk= ik-dS= / ik-sdS , (10) 

Js Js 

where 5* is the aggregate surface, s the unit vector orthogonal to the ag- 
gregate surface, and dS the surface element [see, also, Eq. (|3|]. This cal- 
culation of the collision rate reflects the experimental procedure to measure 
the attachment coefficient: the total mass transfer rate is proportional to 



the integral collision rate. Equation (10) highlights our approximations in 



that the friction coefficient is approximately related to a surface integral of 
the molecular diffusive flux instead of the stress tensor. It is worthwhile 
noting that the linearized Bhatnagar-Gross-Krook equation (BGK), a single 
relaxation-time approximation of the Boltzmann equation, leads to a self- 
diffusion coefficient inversely proportional to the BGK velocity-independent 
collision frequency |21j . 

The proposed methodology suggests that the friction coefficient of fractal- 
like aggregates may be accurately calculated from the solution of the scalar 
Laplace equation for the fluid density (with appropriate boundary conditions) 
without solving the Stokes equations. Once the fluid density has been de- 
termined the collision rates may be easily calculated during post-processing 



via Eq. (10). On the other hand, the Stokes equations (linear equations as 
the Laplace equation) must be solved for the three components of the fluid 
velocity to calculate the drag force. This replacement, in addition to its theo- 
retical implications, has considerable computational advantages that depend 
primarily on the symmetry of the aggregate and the dimensionality of the 
system. Specifically, if axisymmetric aggregates are considered, as in this 
work, the Stokes equations must be solved for the two components of the 
drag force. The Laplace equation, however, must be solved only once for the 
fluid density in two dimensions. 

2.2. Friction coefficient of straight- chain aggregates 

Equations (pi) and (10) will be used to calculate average and individual 



monomer shielding factors and the friction coefficient of straight monomer 
chains via the steady-state collision rate Kk to a single aggregate. We stress 
that Eq. ([s]) is general enough, and easy to implement numerically, to be 
applicable to a generic fractal-like aggregate. A straight-chain aggregate is 
cylindrically symmetric, suggesting that a friction coefficient along the axis 



orthogonal to the axis of symmetry /3^ and parallel to it /3j^ may be defined 

as 



ir 



/3i 



Vk 



/3i 



Vk 



(11) 



The anisotropic collision rates Kj^ are obtained by projecting the diffusive 



kK^ Pi~ " ' kKl 

! obta 
flux Jfc parallel and perpendicular to the symmetry axis, and integrating its 

dp 



absolute value over the aggregate surface, cf. Eq. (|T0|) 

kI^ I 4 dS= f\3k-h\dS= I 
Js Js Js 



Dr. 



dz 



dS 



(12a) 



Ki 



\Jt\ dS 



Hfcl 



-s 1 2 's 1 2 



1/2 



1/2 



t/5 



rf5 



D. 



dp 
dx 



dS 



+ {D, 



dp 
dy 



ql/2 

(1,2b) 



where (in three dimensions and in a Cartesian co-ordinate system) the unit 
vectors are sy = s^, s_l^i = s^, and s_l^i = s^ with the z coordinate along the 
symmetry axis and x,y the co-ordinates perpendicular to it. The absolute 



value is necessary in Eq. (12) to ensure that the collision rate is non-zero. As 



expected, Eq. (11) shows that an isolated spherical monomer has only one 
isotropic friction coefficient, /3^ = jSl 



/3i. 



According to Eq. (11 ) the parallel and perpendicular fluxes to a spherical 
monomer should be calculated. Even though a sphere is isotropic once a 
symmetry axis is randomly chosen the two fluxes differ, K^ ^ K{. Let the 
symmetry axis be the z-axis. Since its choice is arbitrary (due to spherical 
symmetry) 



f\J^\dS= f\Jf\dS= I \Jl\dS 
Js Js Js 



27rDgRip^ . (13) 

'5 JS JS 

However, the perpendicular flux depends on the two fluxes in the two perpen- 
dicular directions. Speciflcally, the (local) molecular fluxes at the monomer 
surface are jj' = Jf = Jicos6 and J^ = (IJ^'P -F |J?p)^/^ = Ji|sin6'| with 
Ji\^.^ji^ = DgPoo/Ri, cf. Eq. (J3|. Then, the analytical evaluation of the 
diffusional fluxes on the surface of the monomer gives 

kI= / jj dS = Dgp^Ri d<p I de\cos9\sme = 27rDgRip^, (14a) 
Js Jo Jo 
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n p27T ptt 

K^ = \J^\ dS = Dgp^Ri del) I d9 sin^ 9 = n^DgRip^ . (14b) 
Js Jo Jo 



Note that, as expected from Eq. (12b), K^ ^ Kf (in fact, K^ > Kf), and 



that Kl is independent of the choice of the symmetry axis, as shown in 



Eq. (13). We stress that even though the parallel and perpendicular diffu- 
sional fluxes as defined differ for a spherical monomer, its diffusion coefficient 
is unique and it does not depend on the choice of the symmetry axis. 

We calculated the friction coefficients (average, and for motion parallel 
and perpendicular to the symmetry axis) of straight chains consisting of up 
to fc = 64 monomers. The finite-element software Comsol Multiphysics [23] 
was used to solve the diffusion equation in cylindrical coordinates (r, z) with 
an absorbing boundary condition on the chain surface. In this co-ordinate 
system j|J = J^ and J^ = J^. We tested the mesh-independence of the 
solutions. The size of the cylindrical computational domain was at least two 
orders of magnitude larger than the corresponding dimension of the chain 
to ensure that the condition p^o = const hold at the computational-domain 
boundaries. The calculated diffusive flux to a single fc-chain, and its com- 
ponents parallel and perpendicular to the symmetry axis, was numerically 
integrated over the aggregate surface to determine the molecule-aggregate 
collision rate. 

Table [T] compares the friction coefficients obtained via the ratio of col- 
lision rates to previous analytical and numerical results. The agreement is 
very good, justifying a posteriori our main approximations. For example, 
explicit calculation of the difference between our results and those reported 
by Dahneke |H] for chains composed oi k = 2, 3, 4, 5, 8 monomers shows that 
the maximum difference is 0.88%, 3.40%, and 1.38% in rjk, rjl, and r/^, re- 
spectively. 

Filippov |6j calculated the friction coefficient of non-overlapping spheres 
via a multipole expansion of the Stokes flow velocity in a series of spherical 
harmonics. The friction coefficient of a straight chain consisting oi k = 8 
monomers was reported to be (3g/(3i = 0.435 (in our notation), in good 
agreement with the calculated 0.428, Table [1} The other values in the second 
column of Table [T] were obtained from reported dynamic shape factors by 
inverting Eq. ([9]). 

Dahneke [8] , based on an interpolation of experimental measurements on 
chains with k < 5, provides extrapolation formulae for the dimensionless drag 
forces felt by a straight chain for motion parallel and perpendicular to the 
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Friction Filippov Happel Brenner Dahneke Collision 
coefficient Ref. [6] Ref. [7J Ref. pj rate 

V2 = N^X 
r/3 = /33//3l 

r/4 = /34//3i 
m = /38//3i 



?ll 



/3|//3i 
/3j//3i 
^i = /3|//3i 
/3S//3i 



0.645 



^3 



4 









0.726 
0.613 
0.550 
0.508 
0.435 



0.716 



0.692 


0.694 


0.569 


0.574 


0.507 


0.507 


0.461 


0.463 


0.390 


0.389 


0.639 


0.633 


0.511 


0.500 


0.442 


0.430 


0.397 


0.385 


0.324 


0.313 


0.719 


0.725 


0.608 


0.612 


0.545 


0.547 


0.504 


0.503 


0.434 


0.428 



Table 1: Comparison of calculated friction coefficients of straight chains of fc = 2 — 8 
monomers (average, and for motion parallel and perpendicular to the chain symmetry 
axis) with previous analytical and numerical results. The values for 77^ — 77^ in the second 
column where obtained from reported values of the dynamic shape factor via Eq. ([9]) . 



symmetry axis. A straight chain composed of k monomers was approximated 
as an ellipsoid of aspect ratio k. In our notation the interpolation formulae 
read 



^k 



(^) 



^||(±)(fc^-l) 
QT^k 



2(A;2-1)±1 



In (k + VF - 1) +5||(x)A; 



(1^) 



where the positive (negative) sign is chosen for the parallel (perpendicu- 
lar) coefficient. The formulae presented in Ref. [S] map to ours via rjl = 
C||,±(0)/(67rA;). We ffited our calculated anisotropic friction coefficients to 
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Eqs. (15). The fits were performed with the Levenberg-Marquardt algorithm 
as implemented in the statistical language R Minpack library [22]. We found 
A^ = 51.07, B± = 0.8637, Ay = 27.45, and Ey = -0.9386 in excellent 
agreement with the fit reported by Dahneke j8] (see, also. Fig. [I]). 

Vainshtein et al. |1] performed an extensive analytical study of the equiv- 
alent mobility radii of oblate and prolate spheroids. Their results for non- 
porous impermeable prolate spheroids with aspect ratio k expressed in our 
notation are 



Vk 



Vk 



Pi 



3k 



F-l 



+ 



/3i 



3k 



2k 



A;2 



+ 



2P 


-3 


(A;2- 


1)3/2 


2A;2 


-1 



\n(k + Vk^-l 



In 



k + VW^i 



-\ -1 



-, -1 



(16a) 



(16b) 



1 {k-'-iy-^r^ \k-Vk^-l, . 

Equations (16), expressed in terms of the mobility radius Eq. ([T]), are iden- 
tical to the equations for the radius of the "equivalent radius" derived in 
Happel & Brenner [7j. 



In the limit of large straight fc-chains Eqs. (15 



16|) show that their per- 
whereas the 



unit-mass friction coefficient tends to zero as Pk — )■ 1/ ln(2fc) 

chain friction coefficient becomes infinitely large as kp^ -^ k/ln{2k). The 

proportionality constant depends on the equations considered and on the 



choice of the constants in Eq. ( 15 ). Moreover, Eqs. ( 16 ) show that the parallel 



and perpendicular friction coefficients tend to /3^ — )■ 2/3]! as A; — )■ oo, whereas 
Eqs. (Q predict /3^ -^ /^I ^±M|| with A^/A\i = 1.86 [see, also, Eq. ^ . 

In Fig. [I] we plot our numerical results (referred to as "Collision rate" , 
circles), the numerical fit based on Eqs. ( jlsj ) ("Collision rate fit", solid line), 
the anisotropic friction coefficients for an impermeable prolate spheroid ac- 
cording to Vainstein et al. [1], Eqs. (16), ("Vainstein et al." , dashed line), 
and the extrapolated results of Dahneke [8] ("Dahneke", dot-dashed line). 
Friction coefficients for motion parallel (left subfigure) and perpendicular to 
the chain symmetry axis are presented. In the right subfigure, squares denote 
results reported by Filippov \Q\ {k = 8) or derived from reported dynamic 
shape factors. As previously remarked the agreement of our calculated fric- 



tion coefficients to those predicted by Eqs. (15) and (16) is very good. 



Happel & Brenner [7] presented a comprehensive analysis of the resistance 
coefficients of two equal-size spheres moving along the line of their centers 
or perpendicular to it. They solved the Stokes equations via the method of 



11 



0.65 
0.55 
^0.45 
$0.35 
0.25 
0.15^ 



(a) 



O Collision rate 

— Collision rate fit 

- - Vainshtein et al. 
■ - Dahneke 




10 20 30 40 50 60 
k 



0.75 
0.65 
^0.55 
H^O.45 
0.35 
0.25 




O Collision rate 

— Collision rate fit 

- - Vainshtein et al. 
■-■ Dahneke 

A Filippov 



10 20 30 40 50 60 
k 



Figure 1: Calculated friction coefficients of straight chains (k — 2 — 64) for motion parallel 
(a) and perpendicular (b) to the chain symmetry axis plotted against the number of 
monomers. Comparison with previous results. 



reflection. Their results for the friction coefficients of two equal-size non- 
rotating spheres as a function of their separation are compared to ours in 
Fig. |2| In particular, for two touching spheres their results /Sg//?! = 0.645 
and ^//3i = 0.716 compare favourably to ours, 0.633 and 0.725, respectively. 
We remark that the comparison is best for touching spheres. 

Table [2] presents calculated dynamic shape factors for the aggregates con- 
sidered in our simulations. The shape factors where calculated via Eq. ^ 
by substituting the appropriate shielding factor iVk^vhvt)- The last column 
presents the results of Ref. [6] determined from a multipole expansion of the 
Stokes velocity field. 

Dahneke [Sj also argued that for a large ensemble of identical straight 
chains undergoing random Brownian rotations the orientation-averaged fric- 
tion coefficient would be 



Pk 



Pi + m 



(17) 



Figure pi compares /3fc//9i evaluated using (^yPi and /3^//3i from the ratio of 
the appropriate collision rates Eq. (12) and Eq. (17) (crosses) to the direct 



calculation of /3fc//3i (diamonds) via Eq. (JSJ) and (10). The agreement is very 
satisfactory. 
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Figure 2: Calculated parallel (a) and perpendicular (b) friction coefficients of two equal- 
sized, non-rotating spheres as a function of their center-of-mass distance z. Comparison 
with the results of Happer & Brenner [7j . 



The mobility radius Rk of straight chains calculated by Eq. ([T]) is shown 
in Fig. |4Fa). As expected, the mobility radii of chains composed of shielded 
monomers are considerably smaller than the mobility radii of ideal chains [for 
which Rk ~ k, Eq. Q]. In addition, their radius of gyration Rg is plotted. 
It is easy to show that the radius of gyration of straight fc-chains is 



1 



Rl = ^Rl{k'-l)+Rl 



(18) 



The last term R\ is an additional term to the usual definition of the radius of 
gyration to ensure that the radius of gyration of a monomer is non-zero. The 



choice shown in Eq. ( 18 ) corresponds to the geometric radius of a sphere [6] ; 
alternatively, the radius of gyration of a sphere (i?i y^3/5) has been used [13] . 
The ratio of the mobility diameter to the radius of gyration is shown in 
Fig. |4Fb). Note that for straight chains the two diameters differ significantly, 
as also noted in Refs. [6l |20]. They argued that the usual assumption Rk ~ 



Rg may be justified for dense fractal-like aggregates with fractal dimension 
greater than 2, but the approximate equality fails for smaller dimensions, in 
agreement with our results shown in Figj4Fb) for straight chains. For large 
straight chains in the limit A; — > oo Eqs. (fTsl) (with the constants determined 



from the fit of the collision-rate friction coefficients) and Eqs. (17, 18) show 
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Number of 


Average dynamic 


Parallel dynamic 


Perpendicular dynamic 


monomers 


shape factor (xk) 


shape factor (xf.) 


shape factor (Xk) 




Collision 


Collision 


Collision 


Filippov 




rate 


rate 


rate 


Ref. [6J 


2 


1.102 


1.005 


1.152 


1.15 


3 


1.194 


1.041 


1.274 


1.276 


4 


1.278 


1.083 


1.379 


1.386 


5 


1.354 


1.126 


1.472 


1.485 


8 


1.556 


1.250 


1.713 




16 


1.975 


1.533 


2.203 




24 


2.312 


1.769 


2.592 




32 


2.601 


1.975 


2.924 




40 


2.862 


2.162 


3.223 




48 


3.100 


2.334 


3.495 




56 


3.322 


2.495 


3.748 




64 


3.531 


2.647 


3.987 





Table 2: Average dynamic shape factor (xk) of straight chains, and their dynamic shape 
factor for motion parallel (Xk) ^-nd. perpendicular {Xk) ^^ their symmetry axis. The last 
column reports results of Ref. [B] . 



that the ratio of the mobility radius to the radius of gyration tends to 



Rk 
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1.05 
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., ln(2A;) 

in good agreement with the slender-body theory prediction (51 [21] 
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V3 
ln(2A;) 



+ 



1 



\n\2k) 



(19) 



(20) 



Equations (16) reproduce the slender-body theory result. 



The connection between collision rates and the average monomer shielding 
factor may be refined to obtain the shielding factor of a monomer within 
an aggregate. The intra-chain (isotropic, per unit monomer mass) friction 
coefficient (3^ of monomer i in a fc-aggregate is defined by 
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Figure 3: Isotropic straight-chain shielding factors 77^ = Pk/Pi calculated via the total 
collision rates from the diffusion simulations (diamonds) and from anisotropic collision 
rates and Eq. (fl7| (triangles) (fc = 2 - 64). 



where Kj^ is the steady-state molecular collision rate of the ith monomer in 
the fc-chain. It is easy to show that 



Vk 



1 '' 
1=1 



Vk 



h = Y.f^ = Y.^h 



(22) 



j=i 



j=i 



.» 



where ft is the total aggregate friction coefficient, fl the total friction coef- 
ficient of the ith monomer, and /i = irtiPi the monomer friction coefficient. 
Equation (21), as Eq. ([s]), can be applied to aggregates of arbitrary shape. 
Table [3] presents position-dependent monomer shielding factors in /c = 5, 8 
straight chains. 



W 









0.597 0.379 0.364 



W 



Vs 



W 



Vs 



Vs 



0.565 0.350 0.325 0.317 



Table 3: Monomer shielding factors in two straight chains, k — 5,8. The superscript 
denotes the relative position of the monomer in the chain, the subscript the total number 
of monomers. 
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Figure 4: Left: Dimensionless straight-chain mobihty diameter Rk/Ri (triangles) and 
radius of gyration Rg/Ri (circles) plotted against number of monomers (fc = 2 — 64). 
Right: Ratio of the mobility radius to the radius of gyration of straight chains as a 
function monomer number {k — 2 ~ 64). 



3. Langevin dynamics of straight-chain aggregates 

Langevin simulations have been used extensively to investigate aggre- 
gate collisional dynamics, see, for example, Refs. [131 ES], and in particular 
the diffusive motion of aggregates. The Brownian motion of an aggregate 
may be described by modelling the Brownian motion of a set of interacting 
monomers held together by strong monomer-monomer interaction forces. A 
coupled set of Langevin equations, each one for a monomer within the aggre- 
gate, is solved to determine the diffusive motion of the whole aggregate. The 
stochastic properties of random force, which describes the effect of molecu- 
lar collisions with the aggregate, are usually determined by assuming that 
the Fluctuation Dissipation Theorem (FDT) hold for each monomer. Ac- 
cordingly, the random force acting on each monomer is determined from its 
friction coefficient. The monomer friction coefficient is usually assumed to 
be independent of aggregate morphology and equal to the average monomer 
friction coefficient /3fc. In particular, in the free draining approximation (i.e., 
for ideal clusters) the monomer friction coefficient is taken to be the fric- 
tion coefficient of an isolated monomer [13]. In the continuum regime, these 
assumptions yield the aggregate diffusion coefficient D^ as given in Eq. (|6|). 
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The previously described determination of the intra-chain shielding fac- 
tor allows the introduction of a monomer-dependent random force. In the 
following we use these intra-chain shielding factors to calculate the diffusion 
coefficient of straight chains from the Langevin dynamics of each monomer 
in the chain. The Langevin equations of motion of the ith monomer in a 
fc-chain is 

mA = F,-/5fmir, + W,(t) , (23) 

where r, is the monomer position, Fj the external force acting on it, and Wj 
the random force that models the effect of molecule-monomer collisions. The 
stochastic properties of the noise are usually taken to be delta-correlated in 
time and space {j,j' = x, y, z; monomers are identified by i, i') 

AJA\\{W^{t)) = 0,{W^{t)W^:{t')) = rAeS,rS{t-t') . (24) 

The monomer-dependent noise strength Fj is determined from the FDT to 
be 

r, = 2/3^*^1 fcsT , (25) 

explicitly depending on the monomer shielding factor. The corresponding 
noise strength in Langevin simulations where the average monomer friction 
coefficient (3k is used becomes F = 2f3kmikBT, i.e., it is independent of the 
monomer position in the chain. 



The force in Eq. (23) models monomer- monomer interactions. It will be 



taken to be conservative 



■Vr,?7, , (26) 



where Ui is the total intermonomer potential the ith monomer feels. The 
monomer potential will be chosen to ensure the stiffness of the straight chain. 
We split potential into a radial part U^'^'^ and an angular part U^^^. The radial 
potential fixes the distance between any two first-neighbor monomers along 
the chain, and the angle-dependent potential U^^^ prevents the chain from 
bending. The radial part is taken to be pairwise additive 

k 

ur" = Y.^^,{n^) ' (27) 

where r^j = |rj — rj|. In our simulations the radial, two-body potential 
Uijifij) was chosen to be the pair- wise additive model potential described in 
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Ref. [13]. The angular potential t/f°^ depends on the angle 0j^j_i^j+i 
by the three consecutive monomers {i — 1, i, i + 1) as follows 



[/f"s = n[l- cos( 



-l,i+l 



vr)] 



k-1 



formed 



(28) 



where f2 is a constant that determines the bending stiffness of the chain. 
Figure [5] illustrates the definition of 0j^j_i_j+i within a k = 5 chain. The 



monomer-monomer potential described by Eqs. (27) and (28) introduces a 



penalty function for any radial or angular deformation of the linear chain to 
ensure that the chain remain straight and rigid during its Brownian motion. 




Figure 5: Definition of the angles used in the angular nionomer-monomer interaction 
potential C/™^ for a linear chain of fc = 5 monomers. 



We performed the numerical simulations of the mobility of straight chains 
with the molecular dynamics package ESPResSo |26J, modified to allow for 
a monomer-dependent random force. The depth of the radial inter-monomer 
potential was chosen deep enough to ensure that monomers remain bound. 
Similarly, the strength of the angular potential Q was chosen to ensure that 
the linear chain did not bend during the simulation {Q = 500/cbT). 

The chain diffusion coefficient Dj^^ was estimated from the time depen- 
dence of the chain centre-of-mass mean-square displacement 



lim {6r, 

t— ^oo 



2 
CM 



it)) = GDt^t 



(29) 



where brackets denote an ensemble average. We followed the motion of the 
aggregate centre-of-mass along 3200 trajectories up to f3it = 100 for chains 
oi k = 5, 8. The chain was placed in the middle of a box of size L with 



L/{2Ri) = 10000. The box size was large enough so that no overall dis- 
placement of the aggregate larger than L/2 was observed. Details of similar 
simulations are described in Ref. flSl. 
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Figure 6: Normalized mean-square displacements of the center-of-mass of straight chains 
consisting of fc = 5,8 monomers {di is the monomer diameter). The linear fits were 
performed for /3ii > 3. 

The dimensionless chain center-of-mass displacements (scaled by the monomer 
diameter di) are shown in Fig. |6]as a function of dimensionless time for the 
two straight chains. The numerical fit of the linear dependence of (Sr^y^) on 
time was performed for (3it > 3. We find D^^/Di = 0.428 and D^^/Di = 
0.319, were Di is the monomer diffusion coefficient. These values are within 
1% of the diffusion coefficient D^/Di = 0.432 and Ds/Di = 0.321 obtained 
from Eq. dGJ) and the shielding factors reported in Table [11 i.e., via the average 
monomer shielding factor calculated from the total collision rates. Since the 
diffusion coefficient is inversely proportional to the friction coefficient these 
results confirm numerically that the aggregate friction coefficient is the sum 



of the fiction coefficient of the shielded monomers, Eq. (22). 

Therefore, the Langevin simulations provide additional support to the 
proposed methodology whereby the diffusion coefficient of straight chains 
may be accurately calculated via the ratio of molecular collision rates. Fur- 
thermore, they suggest that calculations of fc-aggregate diffusion may be per- 
formed by two different, but equivalent, approaches. An average monomer 
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friction coefficient Pk = VkPi niay be used on each monomer in the aggre- 
gate, as in usual Langevin simulations via, e.g., ESPResSo. The FDT is 
assumed to hold for an average monomer and the noise strength is the same 
for all monomers F = 2[5kmikBT. Alternatively, individual monomer fric- 
tion coefficients /3^ {i = 1, . . . fc) are introduced, and the complete set of 
monomer-dependent friction coefficients {/3^ } (or, equivalently, {% }) has 
to be specified. The FDT is, again, assumed to hold for each monomer, but 



the noise strength becomes monomer-dependent Fj = 2/3^ mifcsT, Eq. (23). 
Both descriptions of aggregate Brownian motion are consistent because they 
yield the same diffusion coefficient Dj.. 

4. Conclusions 

We developed a methodology to calculate approximately, albeit accu- 
rately, the friction and diffusion coefficients of generic fractal-like aggregates 
via the ratio of molecule-aggregate collision rates. The collision rates are ob- 
tained from numerical simulations of the molecular diffusion equation with 
an absorbing boundary condition on the aggregate surface. The methodol- 
ogy was validated for straight chains via comparison of friction coefficients, 
isotropic and anisotropic, with analytical calculations and extrapolations of 
experimental measurements. Since the numerical solution of the diffusion 
equation (a scalar Laplace equation) is often more feasible than sophisticated 
treatments of the Stokes equations relatively accurate insights on aggregate 
mobility and monomer shielding factors may be gained from collision-rate 
simulations. 

In addition to the average shielding factor of a monomer in an aggre- 
gate, the collision rate on each monomer was used to define a monomer- 
dependent intra-chain friction coefficient. Langevin-dynamics simulations of 
the diffusive motion of straight chains were performed to argue that within 
our approximation the Fluctuation Dissipation Theorem could be equally 
well applied on the whole aggregate by considering an average monomer or 
on each monomer, as long as the appropriate shielding factors, and conse- 
quently the friction coefficients, are used. The diffusion coefficients of straight 
chains composed of A; = 5, 8 monomers calculated by Langevin simulations 
were found to be in good agreement with the diffusion coefficients calcu- 
lated via the ratio of the molecular collision rates. The Langevin simulations 
depended explicitly on the inter-monomer interaction potential, thereby pro- 
viding a useful technique to investigate the effect of inter-particle forces on 
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agglomeration dynamics and on the shape of the resuhing aggregates. 
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